REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply 
with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To 


15-01-2002 


4.  TITLE  AND  SUBTITLE 


Conference  Proceedinqs 


Irish  Atomistic  Simulators'  Meeting 


17  December  2001  -  18  December  2001 


5a.  CONTRACT  NUMBER 

F61775-02-WF023 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


5d.  PROJECT  NUMBER 


Conference  Committee 


5d.  TASK  NUMBER 


5e.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

National  Microelectronics  Research  Center  (NMRC) 

Lee  Maltings,  Prospect  Row 

Cork 

Ireland 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


EOARD 

PSC  802  BOX  14 
FPO  09499-0014 


11.  SPONSOR/MONITOR’S  REPORT  NUMBER(S) 

CSP  02-5023 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


14.  ABSTRACT 


The  Final  Proceedings  for  Irish  Atomistic  Simulators’  Meeting,  17  December  2001  -  18  December  2001.  The  following 
topics  were  covered  during  this  conference:  computational  materials  science,  electronic  structure  theory,  atomic  scale 
simulation,  multi-scale  modeling,  nanotechnology  simulation 


15.  SUBJECT  TERMS 

EOARD,  Modelling  &  Simulation,  Materials,  Nanotechnology 


20020211  303 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

UNCLAS 


b.  ABSTRACT  c.  THIS  PAGE 

UNCLAS  UNCLAS 


17.  LIMITATION  OF 
ABSTRACT 
UL 


18,  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 
OF  PAGES  Charles  H.  Ward,  Lt  Col,  USAF 
10 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

+44(0)20  7514  3154 

Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39-18 


6th  IRISH  ATOMISTIC  SIMULATION  MEETING 

UNIVERSITY  COLLEGE  CORK 
DECEMBER  17-18,  2001 


We  wish  to  thank  the  following  for  their  contribution  to  the  success  of  this  conference: 

United  States  Air  Force  Office  of  Scientific  Research,  European  Office  of  Aerospace  Re¬ 
search  and  Development. 


CONFERENCE  SCHEDULE 

Venue:  All  talks  will  be  held  in  National  Microelectronics  Research  Centre,  Main  Con¬ 
ference  Room,  Lee  Maltings,  Prospect  Row,  Cork. 

Monday,  December  17 

2:30pm 
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Exchange  coupling  in  CaMn03  and  LaMn03:  UHF 
and  Cl  calculations 


8:00pm  —  Conference  Dinner  at  Proby’s  Bistro,  Proby’s  Quay 

Tuesday,  December  18 

9:00am 
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ABSTRACTS  OF  PRESENTATIONS 


Simulating  the  growth  of  high-k  dielectric  films 
Simon  Elliott 

Materials  Modelling  Group,  National  Microelectronics  Research  Centre,  UCC,  Cork 

Conventional  MOS  transistors  using  Si(>2  gate  oxides  are  still  being  scaled  down,  but  face  a 
severe  restriction  on  minimum  device  size:  they  are  limited  by  the  high  field  gradients  and 
large  tunnelling  currents  across  the  very  thin  gate  layer  of  Si02-  Oxides  such  as  HFO2, 
Zr02  and  A1203  have  high  dielectric  constants  (k)  and  are  being  proposed  to  replace 
Si02,  thereby  promising  comparable  insulating  performance  but  with  thicker  gate  layers. 
We  report  on  a  project  which  aims  to  fill  the  gap  in  understanding  for  the  behaviour  of 
these  new  oxides,  by  multi-scale  simulation  of  the  deposition  process. 


Density  Functional  Theory  study  of  the  absorption  and  diffusion  of  hydrogen 
on  the  [111]  surfaces  of  Ni,  Pd,  and  Pt 

Graeme  Watson 

Department  of  Chemistry,  Trinity  College,  Dublin. 

We  present  gradient  corrected  density  functional  theory  calculations  on  the  adsorption  and 
thermally  activated  diffusion  pathways  of  hydrogen  on  the  (111)  surfaces  of  Ni,  Pd  and  Pt. 
We  find  that  the  variation  of  the  adsorption  energy  as  a  function  of  adsorption  site  shows 
considerable  differences  between  the  three  metals.  For  Ni  and  Pd  the  adsorption  energies 
vary  as  a  function  of  hydrogen  co-ordination,  with  the  3-fold  hollow  sites  the  most  stable 
and  the  1-fold  atop  site  considerably  less  stable.  On  Pt  the  adsorption  energies  for  all  the 
sites  are  similar,  indicating  that  diffusion  across  the  surface  will  be  faster  on  Pt  than  on 
Ni  or  Pd.  The  activation  energies  for  diffusion  have  been  calculated,  with  that  for  Pt  (2-3 
kJ  mol-1)  considerably  smaller  than  for  Ni  or  Pd  (13  kJ  mol-1).  This  is  related  to  the 
experimental  data  on  hydrogenation  activity  and  is  a  possible  contributing  factor  to  the 
differences  observed  in  reactivity. 


Ultrafast  Optical  Excitation  of  Phonons  in  Tellurium 

Sorcha  Healy  and  Stephen  Fahy 

Department  of  Physics,  University  College  Cork 

The  coupled  dynamics  of  photo-excited  carriers  and  the  Ax  and  A2  phonons  are  simulated 
for  ultra-fast  optical  pump-probe  experiments  on  tellurium.  Carrier  diffusion  and  drift 
due  to  electron-phonon  coupling  to  the  A\  phonon  are  considered.  Different  diffusion  rates 
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for  electrons  and  holes  generate  strong  electric  fields  near  the  surface  of  the  photoexcited 
material  on  a  ps  time-scale,  and  these  electric  fields  couple  to  the  infra-red  active  A2 
phonon.  We  examine  the  expected  damping  of  the  Ai  phonon  for  very  large  fluence  of  the 
pump  pulse  and  find  that  it  differs  substanially  from  that  observed  in  experiment. 


Perturbation  theory  for  strongly  frustrated  finite  spin  systems 
N.  P.  Konstantinidis 

Department  of  Physics  and  Centre  for  Scientific  Computation,  Trinity  College  Dublin 
D.  Coffey 

State  University  of  New  York  at  Buffalo,  Amherst,  USA 

Degenerate  perturbation  theory  in  high  orders  is  applied  for  a  series  of  strongly  frustrated 
magnetic  clusters.  An  effective  Hamiltonian  for  the  low  energy  states  is  constructed  as 
a  power  series  in  the  perturbation  parameter  A.  The  series  converge  asymptotically  and 
are  summed  using  Pade  approximants,  providing  the  low  energy  spectrum  and  the  cor¬ 
responding  wavefunctions.  The  first  class  of  systems  investigated  are  ’’toy”  Heisenberg 
models  with  connectivity  similar  to  the  one  of  molecules  of  C,  known  as  fullerenes.  The 
magnetic  properties  of  these  systems  as  functions  of  magnetic  field  depend  strongly  on  the 
connectivity  at  the  classical  level,  and  this  dependence  is  found  to  survive  at  the  quantum 
level.  Furthermore,  the  dependence  of  magnetization  on  the  number  of  sites  is  reflected 
in  the  analytic  structure  of  the  ground  state  energies  and  wavefunctions  as  functions  of 
A.  The  second  system  studied  is  a  model  for  the  magnetic  molecule  V15.  The  model  in¬ 
cludes  a  Heisenberg  term,  and  magnetic  anisotropy  is  included  with  a  symmetry  allowed 
Dzyaloshinskii-Moriya  interaction.  The  properties  of  the  model  are  calculated  as  functions 
of  temperature  and  magnetic  field.  By  comparing  with  data  the  parameters  of  the  DM 
interaction  are  evaluated.  The  magnetic  anisotropy  term  is  found  to  be  an  intrinsic  source 
of  decoherence  for  moments  induced  by  magnetic  fields  even  at  zero  temperature. 


DNA  Triplexes,  a  study  of  their  Hydrogen  bonds 
Martin  Peters 

Department  of  Chemistry,  Trinity  College,  Dublin. 

An  ab  initio  study  of  novel  DNA  triplexes  and  prediction  of  their  hydrogen  bonding  in¬ 
teractions  using  Bader’s  theory  of  ’Atoms  in  Molecules’  is  presented.  We  discuss  how  the 
lapacian  of  the  electron  density  can  characterise  the  triplexes  interactions. 
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First  principles  simulations  of  exchange  coupling  in  La2Cu04:  configuration 
interaction  calculations  and  the  coupling  mechanism 

G.  Zheng  and  C.H.  Patterson 

Department  of  Physics  and  Centre  for  Scientific  Computation,  Trinity  College  Dublin 

The  exchange  constants  of  La2Cu04  have  been  investigated  using  localised  orbital  config¬ 
uration  interaction  (Cl)  calculations  on  clusters.  Exchange  constants  from  ab  initio  cal¬ 
culations  on  clusters  containing  two  and  four  Cu  ions  are  compared  to  estimates  for  these 
quantities  based  on  analysis  of  spin  wave  experiments.  Calculations  show  that  exchange 
interactions  between  the  dx2_y2  orbitals  play  a  very  important  role  in  the  exchange  coup¬ 
ling  of  La2Cu04  and  O  to  da,2_J/2  excitations  have  secondary  importance  in  the  La2Cu04 
ground  state.  The  coupling  mechanism  is  discussed. 


Dynamics  and  free  energy  of  Alumina  surfaces,  semi-empirical  and  DFT  cal¬ 
culations 

Arnaud  Marmier 

Atomistic  Simulation  Group,  Queen’s  University  Belfast 

We  report  an  atomistic  investigation  of  static  and  dynamic  properties  of  the  (0001)  surface 
of  alpha-alumina  (A1203)  at  different  temperatures.  Lattice  dynamics  calculations  are 
performed  within  the  quasi  harmonic  approximation  and  with  two  parameterization  of  the 
shell  model  as  well  as  with  a  DFT  description.  The  simple  models  provide  a  convenient 
benchmark  for  testing  the  convergence  of  different  dynamic  related  properties  with  respect 
to  slab  parameters  (thickness,  vacuum  size)  and  to  vibrational  k-point  sampling.  The 
surface  systematically  undergoes  an  important  relaxation,  whose  details  are  highly  model 
dependent.  The  free  energy  is  calculated  for  different  sampling  of  the  vibrational  Brillouin 
zone.  We  show  that  it  can  already  be  excellently  approximated  with  a  simple  2  k-points 
sampling.  The  analysis  of  the  mean  square  displacement  of  the  atoms  shows  that  the 
surface  atoms  are  between  1.5  and  2  times  more  mobile  than  the  their  bulk  counterparts. 
We  do  not  find  tracks  of  modes  strongly  localized  on  the  surface  A1  atoms  with  the  classical 
model.  On  the  other  hand,  with  the  DFT  model  there  is  clearly  such  a  mode.  However,  it  is 
not  strongly  anharmonic  and  cannot  account  for  the  discrepancies  between  ab  initio  surface 
relaxations  and  experimentally  determined  ones.  This  makes  the  alternative  explanation 
of  surface  impurities  more  likely,  but  so  far  the  experimentalists  have  failed  to  detect  any. 


Exchange  coupling  in  CaMn03  and  LaMn03:  UHF  and  Cl  calculations 
M.  Nicastro  and  C.H.  Patterson 

Department  of  Physics  and  Centre  for  Scientific  Computation,  Trinity  College  Dublin 
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The  equilibrium  structure  and  the  exchange  coupling  mechanism  in  CaMnOa  and  LaMnC>3 
have  been  investigated  using  UHF  calculations  on  the  bulk  crystal  and  Cl  calculations 
on  clusters  representing  the  two  compounds.  In  the  case  of  LaMnC>3  a  series  of  UHF 
calculations  on  a  cubic  idealized  structure  shows  that  independent  orbital  ordering  and 
exchange  terms  contribute  to  the  Hamiltonian.  For  both  systems  Cl  calculations  reveal 
a  strong  dependence  of  exchange  constants  on  the  Madelung  potential  of  the  cluster  and 
on  the  population  of  Mn  d  eg  orbitals.  Analysis  of  the  Cl  wave  functions  shows  which 
exchange  processes  are  important  in  the  exchange  coupling  mechanism.  Predicted  exchange 
constants  are  in  reasonable  agreement  with  experimental  values  estimated  from  the  Neel 
temperature. 


Dielectric  matrix  calculation  and  dielectric  response  eigenfunctions  for  silicon 
and  solid  argon 

S.  Galamic-Mulaomerovic  and  C.  H.  Patterson 

Department  of  Physics  and  Centre  for  Scientific  Computation,  Trinity  College  Dublin 

The  symmetrised  dielectric  matrix  (DM)  has  been  fully  calculated  for  solid  argon  and  sil¬ 
icon  using  Gaussian  orbital  basis  sets.  Calculation  of  the  DM  involves  summation  over 
both  occupied  and  unoccupied  states  and  special  care  has  to  be  taken  while  construct¬ 
ing  the  basis  set.  The  dependence  of  the  DM  elements  on  Gaussian  orbital  bases  will 
be  presented  and  compared  with  published  data.  The  dynamical  dielectric  function  and 
dielectric  response  function  are  calculated  within  the  RPA  for  silicon  and  solid  argon  and 
compared  to  a  full  numerical  calculation  of  these  functions.  The  plasmon  pole-model  for 
the  response  function  is  used  to  obtain  the  plasmon  dispersion  inside  the  first  Brillouin 
zone  and  results  are  compared  with  the  fully  numerical  calculation. 


Energy  optimization  of  correlated  many-electron  wavefunctions 

David  Prendergast  and  Stephen  Fahy 

Department  of  Physics,  University  College  Cork 

In  a  new  iterative  variational  Monte  Carlo  method,  the  total  electronic  energy  is  minimized 
with  respect  to  a  general  expansion  of  the  inhomogeneous  electron  correlation  terms  in  a 
many-electron  wave  function  for  a  periodic  solid.  The  description  of  inhomogeneity  in  such 
systems  represents  a  small  percentage  of  the  correlation  energy,  therefore  determination 
of  the  optimal  variational  wave  function  parameters  is  extremely  delicate.  Our  method, 
based  on  fitting  energy  fluctuations,  is  rapidly  convergent  and  produces  statistically  well- 
determined  parameters  and  correlation  functions,  even  when  a  large  number  of  parameters 
(approx.  3000)  and  electrons  (approx.  200)  are  involved.  Applications  to  the  study  of 
correlation  factors  in  diamond  and  graphite  are  shown. 
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Ab  initio  calculation  of  charged  metallic  surfaces 
Sasha  Lozovoi 

Atomistic  Simulation  Group,  Queen’s  University  Belfast 

Various  properties  of  metallic  surfaces  are  known  to  be  affected  by  the  presence  of  external 
electrostatic  fields  (or  equivalently,  by  appearance  of  surface  charge).  In  addition,  these 
inevitably  arise  in  some  experimental  techniques  such  as  STM,  FIM,  and  electrochemistry. 
We  shall  describe  a  scheme  to  calculate  the  energies  associated  with  charged  surfaces  such 
as  they  could  be  directly  compared.  We  then  demonstrate  how  the  method  works  on 
example  of  the  missing-row  reconstruction  of  Au(110)  and  Pt(110)  surfaces. 


A  molecular  dynamics  simulation  of  an  alkanoate  monolayer  on  water 
Ciara  Bergin 

Department  of  Chemistry,  Trinity  College,  Dublin. 

The  features  associated  with  simulating  an  amphiphile  layer  on  a  water  surface  are  de¬ 
scribed.  A  simulation  of  liquid  water  was  first  achieved  so  as  to  give  satisfactory  values 
of  certain  bulk  properties,  together  with  the  establishment  of  a  liquid-vacuum  interface 
which  disappeared  at  higher  temperatures.  Alkanoate  anions  were  then  sited  with  their 
hydrophilic  regions  (the  polar  COO  group)  below  the  liquid  sin-face,  and  their  alkyl  chains 
directed  into  the  vacuum,  perpendicular  to  the  interface.  The  all-trans  conformations  of 
the  chains  were  approximately  maintained,  but  showed  departures  with  changing  pressure 
and  temperature. 


Ferroelectricity  and  isotope  effects  in  H-bonded  solids:  KDP  and  beyond 
Jorge  Kohanoff 

Atomistic  Simulation  Group,  Queen’s  University  Belfast 

Hydrogen-bonded  molecular  crystals  normally  undergo  an  order-disorder  phase  transition 
at  temperatures  around  100-200  K.  This  is  a  structural  transition  accompanied  by  a  change 
in  the  dielectric  properties  of  the  material,  like  ferro  or  antiferroelectricity.  The  transition 
temperature  is  extremely  sensitive  to  the  isotopic  substitution  of  H  by  D.  This  isotope 
effect  has  been  explained  in  two  ways:  (1)  the  tunneling  model  postulating  that  the  main 
ingredient  is  the  correlated  tunneling  of  the  protons  in  an  effective  double-well  potential, 
and  (2)  the  geometric  effect  which  is  based  on  the  modification  of  the  geometry  of  the 
hydrogen  bonds  upon  isotopic  substitution.  Both  theories  give  results  compatible  with 
experimental  data.  We  have  thoroughly  analysed  the  case  of  KPO4H2  (KDP)  by  means 
of  first-principles  calculations  for  global  and  (correlated)  local  moves  of  the  protons.  Our 
results  indicate  a  predominance  of  the  geometric  effect,  which  is  mainly  connected  to  a 
modification  of  the  lattice  constants.  The  relevant  moves  involve  not  only  protons,  but 
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also  the  displacement  of  heavy  atoms  like  P,  a  fact  that  is  connected  to  the  appearance 
of  ferroelectricity.  The  transition  disappears  if  sufficient  pressure  is  applied.  The  order 
parameter  (polarization),  however,  does  not  vanish  at  Tc,  and  the  protons  remain  off- 
center  in  the  H-bonds  at  a  distance  of  0.1  A  from  the  middle  of  the  bond.  This  result 
appears  to  be  universal  for  the  whole  family.  We  observe  this  feature  in  our  first-principles 
calculations  for  KDP,  and  use  it  to  speculate  about  its  universality. 


An  Exploration  of  Binary  Information  Patterns  in  the  Genetic  Code 
Matthew  Manktelow 

Department  of  Chemistry,  Trinity  College,  Dublin. 

The  universal  genetic  code  employs  the  physical  features  of  nucleotides  to  transmit  inform¬ 
ation.  We  report  the  results  of  qualitative  and  quantitative  analyses  of  these  features  where 
idealised  binary  patterns  are  compared  with  data  derived  from  evolutionary  substitution 
rates.  The  analysis  offers  some  insight  into  the  structure  of  the  genetic  code. 


Surface-Hopping  Dynamics  of  a  Spin-Boson  System 

Donal  Mac  Kernanf,  Giovanni  Ciccotti^’*  and  Raymond  Kapraff’* 

tCECAM,  Ecole  Normale  Superieure  de  Lyon,  Prance 
*INFM  and  Dipartimento  di  Fisica,  Universita  “La  Sapienza”,  Roma,  Italy 
*  Chemical  Physics  Theory  Group,  Department  of  Chemistry,  University  of  Toronto,  Canada 

In  this  talk  the  spin-boson  model  is  solved  within  the  framework  of  quantum-classical 
dynamics  using  our  recently-developed  surface-hopping  scheme.  The  quantum-classical 
equation  of  motion  is  expressed  in  an  adiabatic  basis  and  its  solution  is  constructed  from 
an  ensemble  of  trajectories  which  undergo  non-adiabatic  transitions  and  evolve  coherently 
on  the  adiabatic  surfaces.  Details  of  the  algorithm  for  the  simulation  of  the  dynamics  are 
presented  and  the  method  of  simple  Monte  Carlo  sampling  used  to  evaluate  the  expectation 
values  of  observables  is  discussed.  For  the  spin-boson  model  the  present  implementation 
of  quantum-classical  dynamics  is  exact  and  the  results  of  our  surface-hopping  simulations 
are  in  accord  with  previous  numerically  exact  results  for  this  model. 


Water-mixtures  of  ionic-liquids 
Christof  Hanke 

Atomistic  Simulation  Group,  Queen’s  University  Belfast 

Hydrogen-bonded  molecular  crystals  normally  undergo  an  order-disorder  phase  Ionic  li¬ 
quids  are  basically  molten  salts,  which  are  due  to  the  size  of  the  ions  already  liquid  at 
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temperatures  less  than  100  C.  Some  of  these  liquids  are  vey  hygroscopic  and  change  their 
physical  properties  with  increasing  water  content.  In  this  talk,  some  of  the  properties  of  the 
ionic  liquids  1,3-methylmethylimidazolium  (mmim)  /  Cl  and  mmim  /  PF6  are  discussed, 
both  as  neat  liquids  and  with  different  water-contents. 


Response  function  calculation  for  MgB2 

A.  Rowan,  S.  Galamic-Mulaomerovic  and  C.  H.  Patterson 

Department  of  Physics  and  Centre  for  Scientific  Computation,  Trinity  College  Dublin 

MgB2  was  recently  discovered  to  be  superconducting  up  to  40K.  Its  structure  may  be 
regarded  as  boron  in  layers  with  the  graphite  structure,  which  sandwich  Mg  ions.  Eigen¬ 
functions  of  the  dielectric  response  function  are  used  to  analyse  charge  density  response  to 
an  external  potential  in  MgB2  and  to  identify  the  charge  density  oscillation  mode  likely  to 
be  responsible  for  anomalously  strong  electron-  phonon  coupling  in  this  material. 
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